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We perform direct analysis of mirror mode instabilities from the general dielectric tensor for several model 
distributions, in the longwavelength limit. The growth rate at the instability threshold depends on the derivative 
of the distribution for zero parallel energy. The maximum growth rate is always ~ fcii Vrn and the instability 
is of nonresonant kind. The instability growth rate and its dependence on the propagation angle depend on the 
shape of the ion and electron distribution functions. 

PACS numbers: 



I. INTRODUCTION 



O 

Q_i' Numerous observations of waves in the the Earth magnetosheath, as well as at other planets have stimulated studies of long- 

wavelength and low-frequency modes in high j3 magnetized plasmas. It has been theoretically shown that the features of low- 
c/2 ■ frequency waves in hot plasmas differ significantly from those in cool plasmas, even in the limit corresponding to the usual 
magnetohydrodynamic waves \M- These findings have been subsequently proven by direct comparison with observations [pp. 
However, particular interest to the low-frequency modes in hot plasmas is explained by observations of the mirror modes, which 
£^* were found in planetary magnetosheaths |Q, |], [|, gj , in the solar wind []7|] , in cometary comas [|J g] , and in the wake of Io Q1C| , |ll| ] . 
These modes are nonpropagating zero frequency modes (sometimes considered as the kinetic counterpart of the hydrodynamical 
entropy mode), which are expected to grow in an anisotropic plasma with sufficiently high f3±//3\\ (see, e.g., Hasegawa [113]). 

Usual high amplitudes of observed mirror modes show that they easily achieve the nonlinear regime. At the same time, in 
several cases low-amplitude magnetic field structures with the same properties were observed which may mean that the linear 
and nonlinear mirror mode features are generically related. Yet we do not know so far what makes these modes so ubiquitous 
and what determines their nonlinear amplitudes. 

The early explanation of the mirror instability jl2| ] is based on the simple picture of the adiabatic response of the anisotropic 
Q\ . pressure of magnetized particles. Numerical analyses of the mirror instability in bi-Maxwellian plasmas [ |l3[ [14 , |l5| ] have shown 
that the maximum of the growth rate occurs at k±pi ~ 1 (where pi is the ion thermal gyroradius), which was interpreted as an 
^D , indication on the kinetic nature of the instability. 

At the same time, Southwood and Kivelson [ |lq ] proposed a new explanation of the instability mechanism as a resonant one, 
where the presence of a group of the resonant particles (with v» = 0) plays the destructive role in the mode excitation: the 
growth rate of the instability is claimed to be inversely proportional to the number of the resonant particles. This explanation 
was further reiterated with some modifications by Pantellini and Schwartz [|l7| and Pokhotelov et al. [Q, and used by Kivelson 
and Southwood [n9n for the explanation of the nonlinear saturation mechanism. The analysis of Southwood and Kivelson [[160 
is done in the regime where the phase velocity of the perturbation is much less than the parallel thermal velocity, in other 
words, 7 <C feirUTjii, and therefore, is directly applied only at the very threshold of the instability. At the same time, numerical 
calculations [ |l5| ] show that most important events occur in the range 7 ~ fe||i>Ti||> which is not covered in the previous analytical 
k> studies. 

The previous analytical and numerical considerations of the linear regime of the mirror instability, even in the longwavelength 
limit, are, as a rule, restricted to the usage of the bi-Maxwellian distribution. At the same time particle distributions in collision- 
less plasma ma y su bstantially differ from the Maxwellian. For example, due to the ion heating mechanism at the shock (see, e.g., 
Sckopke et al. pQ]), the magnetosheath ion distributions may well deviate from the bi-Maxwellian. It is therefore of interest to 
study the dependence of the instability on the shape of the ion and electron distributions. 

Yet another argument in favor of the analysis of other distributionsf is that there is no good analytical approximations for the 
dielectric tensor for the Maxwellian plasma in the range |w|/fe|iWri| ~ 1, which forced researchers to consider more convenient 
asymptotics. It is, however, possible to find the shapes of the distribution which allow closed analytical presentation of the 
dielectric tensor in the whole range of phase velocities and make the study of the instability physics more transparent. 

In the present paper we study in detail the dependence of the mirror instability on the shape of the ion and electron distributions, 
using model distribution functions which allow direct explicit analytical calculation of the dielectric tensor. We establish the 



generic relation of the mirror instability with the oscillatory modes when the Landau damping is absent and study the transition 
of damping modes to the unstable regime. We also propose an approximation which is useful for the analytical treatment of the 
instability in the most important range 7 ~ &||^Ti|| m general case. 

The paper is organized as follows. In section II we derive the general dispersion relation in the longwavelength for arbi- 
trary distribution function. In sections DirV| we apply the general analysis to three different distributions. In section M we 



derive the instability condition and the growth rate at the threshold for arbitrary distribution. In section VI we develop a useful 



approximation for the analysis of the bi-Maxwellian-kind distributions in the region of the maximum growth rate. 

II. DISPERSION RELATION IN THE LONGWAVELENGTH LIMIT 

In what follows we will be interested in the longwavelength limit where uj -C fi and kvx <C £1, while maintaining the 
phase velocity finite < ui/k < 00. The last inequality means that the phase velocity does not tend to zero in all propagation 
angle range but it certainly may vanish for particular set of parameters. For simplicity we assume that both ions and electrons 
are Maxwellian in the perpendicular direction, so that (vj_) — 2v 2 ^ ± and (vj_) = 8v^ ± . We also denote (v 2 ) — v 2 -,,, and 

/3|i_L = 2Vyii ± ujp/c 2 £l 2 for each species (subscript i stands for ions and subscript e for electrons). Let us introduce the 
refraction index vector N = kc/w, such that N = (iVj_,0, iV||) = N(sm9,0,cos9). With all this the components of the 
dispersion matrix D,j = N 2 Sij — NiNj — cij take the following form (see Appendix ph: 

Ou = JVf(l-^||-/3x))-l-^-^, (1) 

D12 = 0, (2) 

£> 13 = -JVyM (1 - |(/9|| - /3j.)) , (3) 

D22 = N 2 (l-\ cos 2 0(0|| - ± ) + sin 2 9/3± (4) 

uj 2 - uj 2 
- sin 2 8(rif3 i± Xt + r e Pe±Xe) - 1 - 7^ - jp, 

cj 2 j tan 9 

D 2 3 = ~i— p: {nXi ~ TeXe) , (5) 

D 33 = Ni(l- I (/3|| - f, ± )) - 1 - Sggi ( |1 + *L) (6) 



Till \Pi\\ ^ e ll 



Upi tan 2 9 _ oj 2 e tan 2 9 

r iXi ' 7v> T eXei 



1 n 2 ,w ' n 2 

where /3\\ = /3q + [3 e \\, /3± = (3 l± + (3 e ±, n = ftj_/A||> r e = e ±/0 e \\, and 

r Qf 

X = v n J (u-«u)- g- dv t (7) 

The integration in §ty is taken along the path below the singularity Vu = u. In what follows we shall also assume that cj 2 j /fl 2 ^> 1 
and neglect unity relative to this large parameter (which corresponds to the assumption va <C c, where va — cili/uipi is the 
Alfven velocity). In what follows we also neglect uj 2 /£l 2 — (u) 2 JSl 2 )(rn e /mi). In the above derivation we used cu^/ili — 
— uj 2 /O e in the quasineutral electron proton plasma (this is not correct if any admixture of other charged particles is present). 

In the limit w/O, — > 0(and ui/k finite) the dispersion relation D = dct ||-Dy-|| = splits into two ones. One describes the 
purely transverse Alfven wave (the wave electric field vector in the kBo plane, the wave magnetic field vector perpendicular to 
the external magnetic field) with the dispersion 

u J 2 = k 2 v 2 A cos 2 9{l-\{^-p 1 _)). (8) 

In this wave the absolute value of the magnetic field does not change, but the magnetic field rotates. 
The second dispersion relation reads 



*{Z) = [2 - cos 2 0(/3|, -/3 ± ) + 2 sin 2 80± - 2 sin 2 6{np i±X i + r e Pe±Xe) 
-Z 2 (i A \ cos 2 



Xi Xe 

A|| Pe\\ 



■ 2 a r - - 1 2 n ^ ^ 

sin 9 [nxt - r e Xe\ = 0, 



where we introduced Z = u>/k\\VTi\\ for convenience (ui is complex, in general, so that Z = W + iG), and r, )e = [3i.e±/Pi, e \\- 
Eq. (g) describes elliptically polarized waves with all three components of the wave electric field present, so that in general there 
exists a nonzero component of the wave magnetic field B z = N±E y in the direction of the external magnetic field. These waves 
not only rotate the magnetic field but change its magnitude as well. 

The functions x pl a Y the crucial role in the subsequent analysis. They are defined by the integral containing the distribution 
function f(vn) and cannot be explicitly calculated without particular choice of these distributions. It is common to choose / 
as Maxwellian. In this case x is well-known and tabulated but has good asymptotic expansions only for \Z\ <C 1 or \Z\ 3> 1 
(for electrons Z should be substituted by Z(m e /mi)(vTi\\/vTe\\))- This actually restricts possible analytical considerations of 
the mirror instability only with the range \Z\ <C 1. Yet, numerical analyses show that the most important events occur in the 
vicinity of \Z\ ~ 1 which is unavailable to direct theoretical analysis when Maxwellian is chosen. On the other hand, there 
are vague indications that the qualitative features of long waves (instabilities) in the high (3 more or less sensibly depend on 
the lowest moments of the distribution function (provided it is sufficiently "normal": smooth, no beams, no holes, etc.). It 
therefore makes sense to investigate the dispersion relations for a suitably chosen model distribution so that x can be calculated 
and analyzed in the range \Z\ ~ 1. In what follows we shall use three different distributions for these purposes. The waterbag 
distribution / = Q(vq — vjj)/2zJo will be used for study of the behavior of longwavelength modes and their dependence on 
the plasma parameters in the absence of Landau damping. Here Q(x) — 1 if £ > and Q(x) = if x < 0. The hard-bell 
distribution / = 3(i>q — v 2 )Q(vq — v 2 )/4t>g will allow to include the Landau damping effects, and the Lorentz-like distribution 

/ = (2vq/tt)(vq + v 2 )~ 2 removes the upper limit on the particle velocities. The four distributions (including Maxwellian 
/ = (27ri;^|| ) — x / 2 exp(— v 2 /2v 2 ^,,)) mentioned in this paper are shown in Figure [jj 

III. WATERBAG 

The waterbag distribution / = <d(v 2 — vT)/2vo is somewhat peculiar since the Landau damping is absent. The analysis of 
this distribution allows to establish the generic relation of the instability to nondamping propagating modes. It is easy to find 
that in this case 

1 _ 1 

Xl ~ ~ Z2' Xe ~ 3 - Z 2 fiR (10) 

where (i = m e /mi rj 1/2000, R = (3q//3 e p and v 2 -,* = v 2 /3. In the limit Z = one has d = x{Z = 0) = 1/3. It is 
worth noting that for the Maxwellian distribution d = 1, In this section we use for electrons the approximation of the massless 
bi-Maxwellian (instead of above waterbag, which is used only for ions), for which Xe — 1- The resulting dispersion relation 
(g) is a third order equation with respect to Z 2 with real coefficients. Although this equation can be analyzed directly and even 
solved analytically, graphical representation of the roots is much more convenient. 

Figure [| shows the mode with the highest phase velocity (fast mode) for the case when /3m = (3i±. = j3 e \\ = [3 e ± =0.1 and 
massless bi-Maxwellian electrons. It it worth noting that ions are not isotropic since they are Maxwellian in the perpendicular 
direction and waterbag in the parallel direction. The phase velocity of the fast mode is well above ^Hf-nil so that it has nothing to 
do with the mirror instability. We do not consider this mode in the rest of the paper. We do not consider the Alfven mode either. 
The remaining two low -phase velocity modes are shown in Figure together with lo = k\\Voi (solid line). The upper curve is 
above the resonant region having \u\ > \k\\vn | for all ions. The lower mode is inside the resonant region and would damp if 
there were nonzero df/dvn . 

Figure ^ shows the same two modes but in the case /^m = /3i± = /3 e u — (3 e ±_ — 0.5. In both cases the naive instability 
condition 

K = p ± /0H-2-2/0 ± >O (11) 

is not fulfilled, although in the second case K is closer to the threshold just because of the larger /3j_. There is no much 
difference in the behavior of the two modes for these two cases, except a little stronger decrease of the phase velocities towards 
the perpendicular propagation regime in the higher (3± case. 

Figure & shows the behavior of the two modes in the anisotropic case j3q = /3 e || = 0.1, j3n_ = /3 e j_ = 0.5 (so that K = 1), 
and bi-Maxwellian electrons. The lower mode now remains purely propagating mode for smaller angles (diamonds) but turns 
into an aperiodic instability for larger angles of propagation (stars). The obvious conclusion from Figure Sis that the unstable 
mode has its propagation counterpart for the smaller angles of propagation. The relative growth rate G = T/feiiUTjii ~ 1 is large 
in the whole range of instability, so that the approximation G<1 [n6| is not applicable. 

It is of interest to compare this case with the massless waterbag electrons Xe = 1/3- The corresponding curves in Figure g 
show that there is no instability in this case despite the fact that the threshold (111]) is exceeded. 

Thus, the analysis of the waterbag distribution already shows that (a) there is, in general, the propagating counterpart of the 
mirror instability if Landau damping is absent, (b) the instability threshold and growth rate are sensitive to the details of the 



distribution and not only to the second moment, and (c) the instability is aperiodic, that is, in the unstable range W = and 
G > 0. It can be shown that the last feature is generally valid unless the distribution function is very peculiar (see Appendix y). 

IV. HARD-BELL AND LORENTZ DISTRIBUTIONS 

The waterbag distribution does not allow Landau damping since df/dvn = everywhere. In order to get rid of this restriction 
we consider the hard-bell distribution / = 3(vq — v 2 )Q(v 2 — v?,)/4:Vq, which has nonzero derivative but is is compact (/ = 
for |vi| | > Vq . In this case 



l + ^ m ^-^) 2 + G 



4\/5 (V5 + W) 2 + G 2 

I r r - (12) 

iZ VE-w V5 + W 

H = arctan — h arctan — 

2^5 \ G G 

where Z — W + iG, W and G being real, G > 0, and v 2 . = 5v^,n . The corresponding d = \{Z = 0) = 3/5. The corresponding 

expression for \ e is obtained from (O) by substitution Z — » Z\fpiR. 

In order to analyze non-compact distributions too we shall consider the Lorentz distribution / = (2Vq/7t)(Uq + v 2 ) ~ 2 . In this 
case 

16iZ 3i IZ 3iZ 

* l ~ (i + z 2 ) 3 + T~^z ~ (» - zy + {i-zf 1 ( } 

with v^ = Vyii and d — 3. Again, Xe is obtained by substitution Z — > Z^fjlR. 

We shall also compare the results for these distributions with the bi-Maxwellian. In this case there is no compact analytical 
expression for \ and we use direct numerical calculation. 

In what follows we are interested only in the unstable region. The subparticle mode is expected to be strongly damped in the 
propagation range. The "superparticle" mode is not damped in the hard-bell case and almost not damped in the Lorentz case. 

As the first set of parameters for the unstable regime we choose /3,-ii = /3 e ii = 0.1, /3jj_ = (3 e ± = 0.5, and massless bi- 
Maxwellian electrons \e — 1- Figure^ shows the growth rates for the three distributions. The highest growth rate is for the 
Lorentzian, the lowest is for the waterbag. Figure p\ shows the same growth rates as in Figure R but normalized on kvxi\\ which 
allows to compare growth rates of the modes with the same wavenumber k and different angles of propagation. It is seen that the 
maximum growth rates is achieved approximately at the same angle of propagation rs 60° for all distributions, but the threshold 
angle moves towards more quasiparallel regimes for distributions with stronger tails (Maxwellian and Lorentzian). 

Figure g shows the dependence of the growth rate on fij_ when K = 1 and /3ij_//3 e ii = /3j|i//3 e |i = 1 remain constant. Both 
curves correspond to the waterbag ions and massless bi-Maxwellian electrons. Diamonds stand for the same parameters as in 
Figure |[ crosses correspond to /3%± — 1 and /3m — 0.25. The instability is stronger for higher /3j_, 

In the previous analysis we always used the approximation of massless bi-Maxwellian distribution corresponding to \e = 1. 
Figures [k] and [IT] show the growth rate of the instability when the electron distributions are chosen in the same form as the 
ion distributions. One can see that the waterbag distributions become stable, while the growth rate in the case of Lorentzian 
drastically increases. The ratio of the maximum growth rates shown in Figures |q an d \\A\ roughly corresponds to d e = Xe(Z = 0) 



which shows that the maximum growth rate significantly on electrons (see sections M and VI). 

For other combinations of ion and electron distributions the ratios may be even greater as is seen in Figure [lj, where diamonds 
correspond to waterbag ions and massless bi-Maxwellian electrons, while circles correspond to waterbag ions and Lorentz 
electrons. The (3 parameters are the same for both cases. 

V. NEAR THE THRESHOLD 

It is possible to obtain general results just above the threshold of the instability, where Z = iG — ► +0. For / = f(v 2 ) it is 
easy to find 

f 1 Of 

* = / ~n 7T~ dv \\ 

J iG — V\\ ov\\ 



v \\ df , f df , _ f G df 

v 2 TG^d^ dv ^ = - J de dv " +G J vJT&de 



T-7^7^ dv \\ = - / 77 dv W + G / .,2 I ™ ^ dv \\ (14) 



= - / -Jc dv \\ + nG 7£^n=° = d ~ kG > 



where £ — v 2 /2 is the energy (on the unit mass). Substituting this into (g) and neglecting all terms of the order Z 2 and higher, 
one has 

G = -A/B, (15) 

A = [2 - cos 2 0(/3|| - /3_l) + 2 sin 2 6>/3 ± - 2 sin 2 0(n/3 iX di 

+ r e p e± d e )} ( A + |l ) + sin 2 0( n di - r e rf e ) 2 , (16) 

B = — ^-[2 - cos 2 0(Al - /3j_) + 2 sin 2 6>/? ± - 2 sin 2 9(nP, ± d t 

Pi|| 

+ r-eAade)] + 2sin 2 flnpVL/ti [ -i- + -M - 2sin 2 ^(n* - r e d e ), (17) 

\Pi|| Pe||/ 

where we neglected K e ~ «;, yjm e /m,i. The instability threshold for given is found from the condition (7 = 0, that is, A = 0, 
which gives 

2 + fli - /?|| + sin 2 (9[/3|| + /3j. - 2{nl3 i± di + r e (3 e ±d e ) 

n (18) 

+ (ndi - r e d e ) 2 /(^//3 4 || + de/0e\\)] = 0. 

Since < sin 2 9 < 1, the global instability criterion reads (in the assumption that 2 + /3± > /3|i): 

2( ri (3 l± d t + r e J3 e± d e ) - 2 - 2(3 ± - ( n d t - r e d e ) 2 ( -^ + -^) > 0. (19) 

\Pi|| Pe\\J 

It is instructive to consider several simple cases. We can neglect completely the electron contribution by putting d e — 0, 
which gives the instability criterion in the form 

^>2(1 + A), (20) 



/3|| /3. 



± 



and for the bi-Maxwellian distribution, d = 1, reduces to the naive mirror instability criterion. 
On the other hand, when r e = r^ = /3j_//3|| and d e = di = d, one gets 

^>l + f. (2.) 

P|| P-L 

This condition is harder for more compact distributions (d = 1/3 for waterbag and d = 3/5 for hard-bell) and softer for 
distributions with long tails (d = 1 for Maxwellian and d = 3 for Lorentzian). The global instability condition ( |19{ ) can be 
written in a more symmetric form as follows: 



&1&II (22) 

-2(l +/3j .)(A + A) >0 

VPe|| Pi||/ 

which emphasizes the symmetric role of ions and electrons in the instability onset (cf. Pokhotelov et al. mm). 

Indeed, near the threshold ^/kuVTW <C 1 and the response of both electrons and ions is adiabatic, that is, their inertia does 
not play any role. In these circumstances the mass of the particle is of not importance. Their role in the response to the parallel 
electric field is, however, antisymmetric because of the different signs of the charge: the adiabatic response is obtained from 
eE z — (l/n)(dp/dz) — eE z — ihyp/n = 0. The parallel response plays the crucial role in the instability development. As is 
known the instability occurs because of the breakdown of the local frozen-in condition and efficient drag of particles out of the 
field enhancement into the field depletion region JT6J , O, Hit]. Thus, when the magnetic field is perturbed, B z = Bq + 5B Z , the 
perturbation of the density of the species s is 

Sn s 8B Z (5nl ext) 

— = -fr + ^—i (23) 

where <5n ( s ext) is due to the motion along the field lines. In the adiabatic regime j/kuVx\\ <SC 1 this change can be considered as 
a quasistatic response to the effective potential (j> e g — </> + n s SB z /q s , where <f> is the electrostatic potential, jj, s = (v' 2 _) s /2Bo is 



the average magnetic moment, and q s is the charge of the species. The density response to this effective potential can be found 
from the reduced Vlasov equation 



df s df s d<j> e f f df s 

at " oz oz ovu 



which for d/dt = 7 and d/dz = ik\\ gives 



(5n ( s ext) , f ik\\ df s 



3eff 



dV\\. (25) 



?i 0s 7 j + ik\\v\\ dv\\ 

It is easy to see that in the adiabatic regime near the threshold of the instability, 7^0, this expression reduces to the following 



n 0s Aimosqjrfj 



(26) 



where rp is the Debye length calculated with the parallel distribution function. It is easy to see that r 2 D = v^,u/iu 2 d, where 
d = x(Z = 0). The electrostatic potential <f> can be excluded using the quasineutrality condition Sn e = Siii, which eventually 
gives 



SB, 

Bo 



x T R j +T, 



Ane 2 n (r 2 De + r 2 m ) 



(27) 



where we have taken into account that /1 = Tj_/Bq. Eq. ( |27| ) shows that smaller Debye lengths ro (larger d) result in the stronger 
drag of the particles into the weak field region, that reducing the kinetic pressure response to the magnetic field enhancement 
and supporting instability. Therefore, stronger Debye screening (larger d) would lower the instability threshold, in agreement 
with the found from rigorous calculations. 

From (|15|)-(|l7|) it is easily seen that the growth rate is inversely proportional to «, = —n(df/d£)\ v .. —0, and not to the number 
of particles with vn = (cf. Southwood and Kivelson Jig]). The latter is correct for the bi-Maxwellian distribution since 
(df/d£) oc / in this case. For other distributions this relation may well be wrong. For example, for the waterbag distribution 
(df / d£)\ v ,,—o = and higher order terms should be retained to investigate the behavior near the threshold. It is easy to see 

from ( |14| ) that in this case \ = d ~ ctG 2 , where a = — f v7 2 (df /d£)dv\\ is well-defined. The dispersion relation (g) becomes 
than a first order equation for G 2 , which has one positive solution near the threshold. It is clear that in this case the growth rate 
is determined by the whole distribution and not only by the behavior in vn = 0. 

VI. HYDRODYNAMICAL REGIME 

The previous analysis shows that maximum Z is always of the order of unity or larger, which means that ions no longer 
respond adiabatically to the magnetic field enhancements and their inertia begins to play an important role. This also means that 
it is thermal particles of the ion distribution body with v ~ vrg which are mainly responsible for the instability development 
and not the group of resonant particles with i>|| =0. Figure [k] shows that for some distributions the instability may be very fast 
so that the electron inertia should be taken into account. 

The previous analysis gives a clue to the treatment of the instability in the range of maximum growth rates, where G > 1, 
Let us assume that the distribution function is such that i>|i/(i>||) has a sharp maximum at some v m ~ uj.m. An example of a 
distribution of this kind is the Maxwellian /,; = (l/v27rwj,ii ) exp(— v 2 /2Vj,.ii ) for which there was no good approximation for 
X in the range \Z\ ~ 1 so far. For the aperiodic mirror instability with Z = iG, G > 0, one has 

1 _9f ... f 1 „ df 

'II 



X= -p; ^— dv\\ =- / -=- — 2^11^ — dv \\- (28) 

iG - vn ov\\ " ./ G 2 + v„ " dv» " 



For v m ~ 1 ^ G (v m is normalized on Vt\\) the function (G 2 + v 2 ) l varies slowly in the vicinity of the maximum of 
U11 (df/dvn), so that one may approximate 

1 f Of 1 

* = -^T^r^=^T«;- (29) 

Figure |l^ shows the comparison of the numerically found \ f° r me Maxwellian distribution (v 2 ^ = 2) and Z = iG, G > 
with the approximation d29h. The approximation proves to be very good for G > 1 and is only by the factor 2 smaller at 



G — » 0. Figure [l4| shows similar comparison for Lorentzian. Now the maximum growth rate can be obtained by substituting 
Xi = 1/(G 2 + v^i) in (Q). If G is expected to be high, so that G 2 R/i ~ 1, as it occurs for the Lorentzian e — i distributions in 
Figure 10, the electron inertia should be also taken into account by substituting \ e = l/(G 2 R/i + w^J. If, however, the growth 
rates are relatively modest (as in other cases studied in the present paper), the electrons still respond adiabatically and \ e = d e . 
In the last case (M) turns into a third order equation with respect to G 2 . Finding the maximum growth rate from this equation is a 
technical problem. We shall stop for a while at the physical sense of the above approximation. The dependence of the maximum 
growth rate on v m indicates that the particles with high velocities v ~ Vt\\ are taking part in the process. This is related to 
the dynamic redistribution (closely related to the dynamic Debye screening): if the potential changes quickly the low velocity 
particles do not have enough time to change their position and leave the field enhancements. This redistribution is described by 
the same Eq. ( J25J ) but now j/k\\V<r\\ ~ 1. High velocity particles can leave these regions and reduce the kinetic pressure response 
but their contribution rapidly decreases with the velocity since their number decreases. The increase of redistribution efficiency 
and the decrease of the number of screeners with the velocity increase finds its manifestation in that the main contribution 
belongs to the particles at the maximum of V||/(vi|). Since the redistribution plays the destabilizing role, it can be expected that 
the smaller is v m the higher is the growth rate. This can be seen already from Figure || where the growth rate for Lorentzian 
ions, v 2 — 0.5, is larger than the growth rate for the Maxwellian, v 2 n = 2 (with the same massless Maxwellian electrons). 
Figure O^ shows the comparison of the growth rates obtained with the proposed approximation for several v 2 ^ — 2 (diamonds), 
1.5 (crosses), 1 (triangles), 0.5 (circles), and massless Maxwellian electrons. The parameters chosen are (3i± = /3 e ± = 0.5, 
An = /3 e || ■ As expected the decrease of v m results in the increase of the maximum growth rate. 

Finally, Figure [Fg shows the comparison of the growth rates obtained directly and with the above approximation for 
Maxwellian (diamonds and crosses) and Lorentzian (triangles and circles), for the same parameter set. The agreement is quite 
satisfactory. 

VII. CONCLUSIONS 

We have derived the most general dispersion relation for longwavelength modes in hot plasmas. We have derived the general 
mirror instability condition for arbitrary ion and electron distributions and growth rate of the instability near the threshold. 
The instability threshold depends not only on the plasma species (3 but also on another integral characteristic of the distribution 
function d = j v7 l {d f / dv\\)dv\\ for both species. Larger d corresponds to smaller Debye length. Smaller Debye length, in turn, 
corresponds to stronger response of the density to the perturbations of the potential, which allows stronger density depletions 
in the regions of the magnetic field enhancements. Therefore, the kinetic pressure response to the magnetic pressure buildup 
weakens. Hence, the larger is d the lower is instability threshold. The near-the-threshold growth rate is inversely proportional to 
df/d£, where £ — v 2 /2 is the parallel energy. 

The mirror instability is always aperiodic and { / y/k\\Vxi\\)m3x ~ 1 (and sometimes substantially greater). Maximum growth 
rates are normally determined by v m i such that vndfi/dvn has a sharp maximum in v» = v m i, and d e (if the instability is very 
strong v me takes the place of d e ). This is related to the dynamic redistribution in which the thermal particles participate. Growth 
rates are higher for distributions with tails and lower for compact distributions (those, for which / = if \v\\ \ > vq, where vo 
is some upper limit). For noncompact distributions the maximum growth rate is larger for smaller v m , which corresponds to 
the weaker dynamic screening of the parallel electric field. For the distributions analyzed in this paper the behavior of d and 
v m correlates (d increases when v m decreases) since all these are single-parameter distributions. For more general distributions 
the behavior of d e and v m may be uncorrected. It is also worth noting that it not, in general, any specific group of particles 
which are responsible for the instability development. Compare, for example, two similar distributions (velocity normalized on 
the thermal velocity v T \\): /i = (2/tt)(1 + v 2 y 2 with d = 3 and v 2 n = 0.5, and f 2 = (V2/tt)(1 + w|) _1 with d = 1 and 

v 2 n — 1, While the behavior of the two is similar for v\\ — and v\\ — > oo (the only difference is the factor v2), the first one is 
expected to be more unstable because of the three times stronger Debye screening. At the same time the behavior of the second 
distribution near the threshold should be close to that of the Maxwellian, d = 1, despite the very different suprathermal tails and 
(df/d£)\ vll =o- 

We have also proposed a useful approximation for the dielectric function in the range G/k\\Vj-i\\ ^ 1 for distributions with 
sharp maxima of uii (df/dvu ) (Maxwellian as one of such distributions). This approximation proves to be quite satisfactory for 
Maxwellian type distributions and allows to study analytically the instability behavior in the maximum growth rate range. 
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APPENDIX A: GENERAL EXPRESSIONS 



We start with the general expression for the dielectric tensor in the following form: 



where the summation is on the species and 
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The expression for r]ij is well-known (see, e.g., Hasegawa JI4]) 
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Here J n — J n (x), x = k±p = k±v±/fl, and J' n = dJ n /dx. 

For the analysis in the low-frequency range lj/H <C 1 let us write 
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Let also f = fi(v±)f2(v 2 ), and denote (...)= /(. . . )fdvj, where j =_L, 
One has 
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Now, up to 51 2 one obtains 
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Using in Eqs.(A8)-(A13) the following summation rules 
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one obtains eventually the following general expression for A^ in the limit of u>, k\\V\\ <C ft when expanded up to the second 
order in £/ft: 
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It is possible to get rid of the series in ( A25 )-( A27 ) using [Ell] 
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(A29) 

(A30) 



This may be useful for calculations in the regime k±v±/fl ~ 1. 

The general dispersion relation is obtained from the determinant det \D\ = 0, where Dy — N 2 Sij — NiNj — tij, that is, 
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The polarization should be found from the equations 
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Eq. ( A.37 ) provides the ratio of the electric field components . In order to translate that into the magnetic polarization one has to 
use the relation B = k x E/w. In order to find the density perturbations one has to use the current conservation as follows 



5p = k-j/w, 



where 



so that one has eventually 
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Further simplifications are possible in the longwavelength limit. 



APPENDIX B: LONGWAVELENGTH APPROXIMATION 



In this appendix we provide general expressions for the dielectric tensor in the longwavelength limit k±v±/fl <C 1, where 
J±i = ±k±v±/2£l, J = 1 — k 2 _v 2 _/2n, 2 , and higher order Bessel functions may be neglected. In this limit one has 
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2/„,2 \ /^2 2 



where N = fec/w, iVj. = k ± c/uj = Nshx6, N\\ = ^yc/w = iVcostf, /3|, = 2uj 2 p {v 2 ) /c 2 Vt 2 , 0j_ = 0J 2 p {v\) /c 2 Vt 2 , and 



denotes usual averaging over the distribution. Here also £ = u — v\\, where u — w/feii. The last term in (B6) is given for 



completeness. In the limit used in this paper, lu/Q — > and uj/k finite, it should be neglected. Throughout the paper we also 
assume to vi ^> $7j. 

APPENDIX C: APERIODIC NATURE OF THE MIRROR INSTABILITY 

In order to show that the mirror instability is aperiodic we analyze the behavior of the roots of (^) when the (3 parameters 
are changed. In the waterbag case the transition from the stable to the unstable regimes occurs when W = 0, G = and for 
the mode whose phase velocity is less than the highest particle velocity, RcZ < f|| max ("subparticle" mode), that is, in the 
resonant region. In the general case, where Landau damping is nonzero, in the resonant range every propagating wave having 
W y^ has also nonzero damping rate G < (we assume that there are no other kinetic instabilities in the mirror-stable 
region). By continuously changing the plasma parameters (e.g., the anisotropy ratio (3±/f3\\) we can bring the system into the 
unstable regime. Assuming continuous dependence of W and G on the plasma parameters we see that it is impossible for the 
"subparticle" mode with W ^ to transform into the unstable mode, since G < and cannot be made positive continuously. 
Thus, the only way to do that is to go through W = 0, G = 0. 

Let us now consider the vicinity of the transition to the instability, \Z\ <C 1. In the most general way, expanding \ m powers 
of Z one gets: 

X = 

j /J — v\\ av± 
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provided (df/dv\\)\ v ,. = o =£ 0. The quantities d and k are defined in (14). It is easy to see that (|9|) is the first order equation 
for iZ (with real coefficients) in the lowest order on \Z\ <C 1, which means that there is a simple (one and only one) aperiodic 
root in the vicinity of Z = 0. Such aperiodic solutions cannot be converted into non-aperiodic ones by continuous change of the 
plasma parameters, for the same reason as above. Therefore, the unstable solutions must be aperiodic. 

The function ^f(Z), defined in (g), is an analytical function of Z = W + iG and a continuous function of its parameters 
(3 and 9. Let us consider how Z moves from the lower half-plane (stable regime) to the upper half-plane (unstable regime) 
with the change of (3 and 6 = const. The transition to instability occurs, in general, in the vicinity of Z = where &(Z) = 
*(0) + (d^/dZ)\ z= oZ = A + BZ (see sec. |v|). In the transition point A = 0. Using ([B])-([l7|) it is easy to show that in the 
transition point B > (provided k > 0, this condition being violated if df /d£ > at v\\ = 0, corresponding to the regime 
of two-hump instability), so that in the vicinity of the transition point A > corresponds to the stable regime, while A < 
corresponds to the instability. Because of the continuity, in the whole instability range A < 0. 

Let us show now that (ffl always has a solution Z = iG, G > in the unstable range. Indeed, ^(0) < as is shown above. 
On the other hand, if G — > oo one has x ~^ 1/G 2 and ^(oo) > 0. This means that there exists G > such that ^(G) = 0. 

In the absence of kinetic instabilities, in the stable regime all roots of (ffl with nonzero W are either in the lower half-plane 
(Landau damping or nonpropagation) or at the real axis (if (d f / dvn)\ v ,,—w — 0). In the first case no root can cross the real axis 
except at W = 0, when the parameters are changed continuously to bring the system in the unstable regime. As can be seen 
from (p"5|)-(|T7|) there is only one root crossing the real axis at this point, provided df/d£ < 0. Therefore, there is only one root 
in the upper half-plane and it is purely imaginary. 

If df/d£ = there are two or more roots in the vicinity of Z — (depending on the behavior of /) but only one is positive, 
G > 0. Since in this case the analytical continuation through Z = into the lower half-plane is straightforward (no pole at 
I'll = 0) other roots correspond to damping solutions, and there is again only one root in the upper half-plane. 

Finally, let us consider the case where there are roots with G = and W — Wo ^ 0. Such situation can occur when 
(df/dvn ) = in isolated points or in an interval (as for the compact waterbag and hard-bell). In the first case the imaginary part 
of Z is negative for W close to Wo, so that the continuous change of parameters does not bring the root to the upper half -plane. 

In the second case the continuous change of parameters leaves the root on the real axis until it enters the range where 
(df/dv^^OoiW = 0. 
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FIG. 1 : Waterbag (solid line), hard bell (dashed), Lorentz (dotted), and Maxwellian (dash-dotted) distributions. 
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FIG. 2: Phase velocity of the fast mode as a function of propagation angle for the case of the waterbag distribution with /3jm = /3,j_ = /3 e 
/3 e ± =0.1 and massless bi-Maxwellian electrons. 
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FIG. 3: Phase velocity (diamonds) of the two low- velocity modes as a function of propagation angle for the case of the waterbag distribution 
with /3jii = /3i± = /3 e ii = f3 e ± — 0.1 and massless bi-Maxwellian electrons. The solid line is w = y/alkuv Ti \\ = fciiVo;- 
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FIG. 4: Phase velocity (diamonds) of the two low- velocity modes as a function of propagation angle for the case of the waterbag distribution 
with ,0j|| = /3i± = /3 e ii = f3 f ,± — 0.5 and massless bi-Maxwellian electrons. The solid line is u> = ^/aikuVTiU = k\\Voi- 
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FIG. 5: Behavior of the two low- velocity modes as a function of propagation angle for the case of the waterbag distribution with /3i\\ = /3 e \\ — 
0.1, /3i± — /3 e x = 0.5 (so that K — 1), and massless bi-Maxwellian electrons. Diamonds mark the modes in the range where their frequencies 
are purely real, stars show the growth rate of the aperiodic instability. 
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FIG. 6: Behavior of the two low- velocity modes as a function of propagation angle for the case of the waterbag distribution with /3i\\ = j3 e 
0.1, fii±_ = f3 e ± = 0.5 (so that K = 1), and waterbag electrons, d e = 1/3. There is no instability. 
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FIG. 7: Growth rates for the mirror instability in the case of f5i\\ = /? e || = 0.1, /3i± = /3 e ± = 0.5, and massless bi-Maxwellian electrons 
d e = 1, and four different distributions: waterbag (diamonds), hard-bell (crosses), Lorentz (circles), and bi-Maxwellian (triangles). 
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FIG. 8: Growth rates for the mirror instability in the case of /3i\\ = j3 e \\ = 0.1, /3i± = /3 e ± = 0.5, and massless bi-Maxwellian electrons 
d e = 1, and four different distributions: waterbag (diamonds), hard-bell (crosses), Lorentz (circles), and bi-Maxwellian (triangles), normalized 
onkv Ti \\. 
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FIG. 9: Dependence of the growth rate on /3± with K — 1 for waterbag ions and massless bi-Maxwellian electrons: diamonds correspond to 
A|| — Pe\\ = 0.1, pix — /3 e x = 0.5, crosses correspond to (3^ = /3 e ii = 0.25, /3i± = /3 e ± = 1. 
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FIG. 10: Growth rates for the mirror instability in the case of ft || ~ P e \\ — 0.1, /3i_i_ = /3 e J_ = 0.5, and three different combinations: hard-bell 
ions and electrons (crosses), Lorentz ions and electrons (circles), and bi-Maxwellian ions and (massive) electrons (triangles). 
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FIG. 1 1 : Same as in Figure ll(j but normalized on kv 
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FIG. 12: Growth rates for the mirror instability in the case of fii\\ = /3 e \\ = 0.1, fin. — fie± — 0.5, waterbag ions and two different electron 
distributions: massive bi-Maxwellian (diamonds) and Lorentz (circles). 
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FIG. 13: Approximation of x(G) for the Maxwellian distribution. The numerically calculated x(G) (solid line) is compared to \ — (G 2 +2)~ 
(crosses). 
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FIG. 14: Approximation of x(G) for the Lorentzian distribution. The numerically calculated x(G) (solid line) is compared to x = (G 2 + 
0.5)~ (crosses). 
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FIG. 15: Growth rates for the mirror instability in the case of /3i\\ — /3 e || = 0.1, (3i± = /3 e ± = 0.5, calculated with the approximation 
X% = V(C 2 + Vm), for several v^ — 2 (diamonds), 1.5 (crosses), 1 (triangles), and 0.5 (circles). The electrons are massless Maxwellian. 
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FIG. 16: Comparison of the growth rates for the mirror instability in the case of /3jii = /3 e ii = 0.1, /3;x = /3 e x = 0.5, calculated directly 
and with the approximation Xi = 1/(G 2 + uj,), for Maxwellian (diamonds and crosses, respectively) and Lorentzian (triangles and circles, 
respectively). The electrons are massless Maxwellian. 



